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Abstract 



We discuss possibilities to investigate the effects of CP (and T) violation 
in the lepton sector in neutrino oscillation experiments. We consider the 



' effects of CP violation in the framework of two schemes of mixing of four 

^ massive neutrinos that can accommodate the results of all neutrino oscillation 

experiments. Using the constraints on the mixing parameters that follow from 
the results of short-baseline neutrino oscillation experiments, we derive rather 
strong upper bounds on the effects of CP violation in ^v^^^^vl transitions in 
long-baseline neutrino oscillation experiments. We show that the effects of 
CP violation in 'z/^^V^ transitions in long-baseline oscillation experiments 
can be as large as is allowed by the unitarity of the mixing matrix. The 
matter effects, which complicate the problem of searching for CP violation 
in long-baseline experiments, are discussed in detail. We consider the T-odd 
asymmetries whose measurement could allow to reveal T and CP violation in 
the lepton sector independently from matter effects. 
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I. INTRODUCTION 



The violation of CP invariance is one of the most important problems in particle physics. 
So far, CP violation has only been observed in the K^K^ system Many future exper- 
iments will investigate the effects of CP violation in decays of B mesons and are aimed to 
reveal the origin of CP violation in the quark sector (for recent reviews see, e.g., Ref. 0). 
In the Standard Model of electroweak interactions CP violation resides in a phase in the 
Cabibbo-Kobayashi-Maskawa mixing matrix 0]. In the lepton sector an analogous mixing 
matrix is expected to exist if neutrinos are massive and, consequently, it is plausible to 
assume the presence of CP-violating phases in the neutrino mixing matrix as well. 

At present there are some indications that neutrinos are massive and mixed particles 
coming from the results of solar neutrino experiments (Homestake |Q, Kamiokande [Q, 
GALLEX ip, SAGE and Super-Kamiokande |^), of atmospheric neutrino experiments 
(Kamiokande d, 1MB Soudan 



and Super-Kamiokande |T^) and of the accelerator 
LSND experiment [|T^. The analysis of the data of these experiments in terms of neu- 
trino oscillations indicate the existence of three different scales of neutrino mass-squared 
differences: 



Am^ ~ 10" 



eV^ [14] or Amf^^ ~ IQ-^OeV^ [0 



5 X 10-3 eV^ PI 



0.3< Am^g^D <2.2eV2 [0. 



:i.2) 

:i.3) 



The two possibilities for Am^^^ correspond, respectively, to the MSW [|T^] and to the vacuum 
oscillation (see Refs. |17- [l9|] ) solutions of the solar neutrino problem. At present there is no 
information on CP violation in the lepton sector. 

Here we consider possibilities to reveal effects of CP violation in neutrino oscillations 
in schemes of neutrino mixing that can provide three independent neutrino mass-squared 
differences. These schemes are based on the assumption that the flavour neutrino fields are 
superpositions of four massive neutrino fields. This means that the neutrino mixing matrix 
can contain CP-violating phases and effects of CP violation in the lepton sector could be 
observed in neutrino oscillation experiments. 

In principle, short-baseline (SBL) accelerator neutrino oscillation experiments could be 
important sources of information on CP- violation in the lepton sector, but in the case of four 
massive neutrinos only the largest mass-squared difference Am^ = AmLgND is relevant for 



SBL oscillations and the effects of CP violation cannot be revealed in SBL experiments [20 



This fact is discussed in Section |T|, where we review some general aspects of CP violation 
in neutrino oscillations. In the rest of the paper we discuss the effects of CP violation that 
can be expected in future accelerator long-baseline (LBL) neutrino oscillation experiments 
(K2K 



g, MINOS m, ICARUS [H and others [||]). 
In Refs. [p5| , p6i we have shown that among all the possible schemes with four massive 
neutrinos only two can accommodate the results of all neutrino oscillation experiments (see 
also Ref. [^). These two schemes are presented in Section |T|. In Section |V| we apply to 
these schemes the general methods presented in the Appendices that allow to obtain 
limits on the parameters that characterize the CP-odd asymmetries in different LBL channels 
from the exclusion plots obtained in SBL experiments. We show that in the schemes under 
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consideration there are rather severe constraints on the parameter /g^ that characterizes CP 

violations in the ''v^^^vl channels in vacuum, whereas the parameter I^^. that characterizes 

CP violations in 'z/^t^'z/^ LBL transitions in vacuum can reach its maximum value determined 
by the unitarity of the mixing matrix. 

The possible effects of CP violation in LBL neutrino oscillation experiments have been 
discussed recently in the literature |28-^] in the framework of three- neutrino schemes that 
could accommodate the results of some, but not all, neutrino oscillation experiments (in 
particular, in these schemes the solar neutrino problem cannot be explained with neutrino 
oscillations). In Section we apply the methods presented in this paper in order to obtain 
the limits on the parameters that characterize the CP-odd asymmetries in different LBL 
channels in the framework of these schemes of mixing of three massive neutrinos. 

In Section |V| we discuss the implications of matter effects for the possibility to observe 
CP violation in LBL experiments. Matter effects can be large and they represent a serious 
problem for the investigation of CP violation in LBL experiments because the interaction 
of neutrinos and antineutrinos with matter is not CP-symmetric. In order to extract the 
CP-violating phases of the mixing matrix from the measured asymmetries it is necessary to 
have detailed information on the absolute values of the elements of the mixing matrix and 
on the neutrino mass-squared differences. To avoid this problem, we consider the T-odd 
asymmetries whose measurement could reveal CP violation in the lepton sector indepen- 
dently from the presence of matter (the matter contribution to the effective Hamiltonian is 
T-symmetric) [pl|. Measurements of such asymmetries may be possible in the future, for 



example, with neutrino beams from muon colliders |32,33 



II. NEUTRINO OSCILLATIONS AND CP VIOLATION 

In accordance with the neutrino mixing hypothesis (see, for example, Refs. p!7|-[T9|), a 
left-handed neutrino field VaL is a mixture of the left-handed components vul of the (Dirac 
or Majorana) fields of neutrinos with definite masses m^: 

VaL = '^UakVkL with « = 6, /i, T, S, . . . , (2.1) 

k 

where U is the unitary mixing matrix. Here > 3 and z/^^, . . . are possible sterile neutrino 
fields. The mixing in Eq. (|2.1|) implies that the transition probabilities in vacuum of neutrinos 
and antineutrinos with momentum p at a distance L of the neutrino detector from the 
neutrino source are given by 

2 

(2.2) 

2 

(2.3) 

where Am^^ = m\ — m\ (we take mi < m2 < . . .). From Eqs.( p.2|) and (|2.3| ) it follows that 
the transition probabilities of neutrinos and antineutrinos are connected by the relation 



p 



lik U^j, exp 



2p 



P- 



*U0 



Uak exp 



2p 



3 



Pua^Ufj — Pup^Ua ■ (2-4) 

This relation reflects CPT invariance. 

If CP invariance in the lepton sector holds, then there are conventions for the arbitrary 
phases such that in the case of Dirac neutrinos we have 

Uak = K,, (2.5) 

whereas in the case of Majorana neutrinos we have 

U^k = -U:,Vk, (2.6) 

where rj^ = ±i is the CP parityQ of the Majorana neutrino with mass (see, for example, 
Ref. |18[)- It is obvious that the CP parities r]k do not enter in the expressions for the 
transitions amplitudes. Hence, in both the Dirac and Majorana cases, CP invariance implies 
that [0 

Let us introduce the CP-odd asymmetries 

Da;l3 = Pu^^vp - Pv^^up ■ (2.8) 

From CPT invariance it follows that 

Da;l3 = -Dp,^ . (2.9) 

Furthermore, from the unitarity of the mixing matrix we have 

Y,Da;P = Q. (2.10) 

We observe that in the case of transitions among three flavour states (i^e, t'/i, Vt) the CP 



asymmetries satisfy the relations [20 



D,,^ = D^,r = Dr,e, (2.11) 

which follow from Eqs.( |2.9| ) and ( |2.10| ). 

In the general case of mixing of an arbitrary number of massive neutrinos, the transition 
probabilities are given by 

P(_) ,_,=Y^\U^^\^\Up^\' + 2Y,MUo.iU;^KkUpk]cos^^±-D^,p, (2.12) 



^ The CP parities of Majorana neutrinos could be important for neutrinoless double beta-decay; for 
example, if the f^'s have different CP parities, their contributions to the amplitude of neutrinoless 
double-beta decay could cancel each other [34|. 
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where the plus (minus) sign apphes to neutrinos (antineutrinos). The expression for the 
asymmetries is 



Da;(3 = Iaf3;jk Sm , (2. 13) 

k>j P 

with 

= 4 Im [f/„, U;^ U:, U^k] . (2. 14) 

These parameters are invariant under rephasing of the neutrino mixing matrix and (apart 
for a factor 4) are the analogues in the lepton sector of the well-known rephasing-invariant 
parameters in the quark sector |]36|- 0^] . In Sections ^VHV^ we will derive the constraints on 



the parameters lap-jk which follow from the results of neutrino oscillation experiments in 
the framework of schemes with four and three massive neutrinos. 

CP violation in the lepton sector can be observed in neutrino oscillation experiments only 
if at least one of the terms of the sum in Eq. ( |2.13| ) does not vanish because of the averaging 
over the neutrino energy spectrum and the size of the neutrino source and detector. 

If only one mass-squared difference (denoted by Am^) is relevant for short-baseline neu- 
trino oscillations, the neutrinos can be divided in two groups z/i, . . . , z/^ and z/^+i, with 
masses rrii < . . . < rrir and m^+i < . . . < m„, respectively, such that in SBL experiments 

<^ 1 for j,k<r or j,k>r, (2.15) 



2p 



whereas 



Amij ~ Am^ for k > r and j < r . (2.16) 
In this case, for the CP asymmetries in SBL neutrino oscillation experiments we obtain 

^1T'= (E + E + E) W=in^- E W=in^ = 0. (2.17) 

\r>k>j k>j>r k>r>j/ k>r>j 



The last step follows from the definition ( p.l4|) and the unitarity of the mixing matrix. 



Consequently, it is necessary to consider neutrino oscillations in LBL experiments in order 
to have some possibility to observe effects of CP violation in the lepton sector. 

In the rest of this paper we will consider schemes with four and three massive neutrinos, 
in which only the largest mass-squared difference Am^ is relevant for neutrino oscillations 
in SBL experiments, having a value in the wide range 

IQ-^eV^ < Am^ < 10^eV^ (2.18) 

which include the interval ( p..3|) allowed by the results of the LSND experiment. 
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III. FOUR MASSIVE NEUTRINOS 



All existing indications in favour of neutrino oscillations can be accommodated by a 
scheme with mixing of four massive neutrinos [^ , ^ ,p5|-p7| . In Refs. p5| , p6| we have shown 
that from the six possible spectral schemes of four massive neutrinos, which correspond 
to three different scales of mass-squared differences Am^^-, only two schemes are compatible 

In these two schemes the four neutrino 



with the results of all experiments (see also Ref. |27 
masses are divided in two pairs of close masses separated by a gap of 



leV: 



atm solar 

(A) mi < m2 <^ 7723 < mi" 

^ V ' 

LSND 



solar atm 

and (B) mi < m2 <^ m-t < ml 
^ V ' 

LSND 



(3.1) 



In scheme A, Am^^ is relevant for the explanation of the atmospheric neutrino anomaly and 
Am43 is relevant for the suppression of solar z/g's. In scheme B, the roles of Am^i and Am43 
are reversed. 

Let us define the quantities (with a = e, fi, r, s) as 



5^ \Uak\ 



(3.2) 



fc=l,2 



Taking into account the results of SBL neutrino oscillation experiments and those of solar 
and atmospheric neutrino experiments, in the two schemes A and B the parameters Cg and 



are constrained by p5| , p6 



(A) 
(B) 



Ce < di 



Ce > 1 



> 1 



where 



«° = ^ (l - V^l^^) (« = e,/x) 
and S^.Q, is the experimental upper bound for the oscillation amplitude 



Boi\a 4Cq,(1 Co 



(3.3) 
(3.4) 



(3.5) 



(3.6) 



in SBL disappearance experiments. The values of Og and a° obtained, respectively, from 
the 90% exclusion plots of the Bugey |^ Ve reactor experiment and of the CDHS |^ 
and CCFR [43| v^,^ Vy, accelerator experiments are given in Fig. 1 of Ref. [0]. From that 



figure one can see that is small (a^ < 4 x 10^^) for Am^ in the wide range ( p.l8| ) and a° 
is small (a^ < 10"^) for Am^ > 0.5 eV^. In the following we will use also the upper bounds 
y4°.^ for a = /i and /5 = e, r on the SBL oscillation amphtudes p6| 



A 



E 

fc=l,2 



E 

fc=3,4 



(3.7) 
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which are obtained from 90% CL exclusion plots of the BNL E734 [||], BNL E776 |6 



and CCFR m 



(-) (-) 



(-) 



(-) 



Up appearance experiments and of the FNAL E531 |^ and CCFR 

experiments. The result of the LSND |T^ — > z/e appearance 
experiment, which is crucial for the arguments in favour of the 4-neutrino schemes A and 
B p3|^7[, will be of importance in the further discussion mainly through the allowed range 
(|OD of Am2. 

It can be seen from the neutrino mass spectra (|3.] 



that the expressions of the oscillation 
probabilities in scheme B follow from the corresponding expressions in scheme A through 
the exchange of indices 



1, 2 t^3, 4. 



(3.8) 



Since this permutation of indices transforms the conditions ( p.3| ) into 1 



Ce, < 



and 



> 1- 



which are identical with those of Eg. ( |3.4|) , it follows that the schemes A and 



B are indistinguishable with neutrino oscillations Note, however, that these schemes 
could in principle be distinguished, for instance, in (/5/3)oi/ decay or with the measurement 
of the high-energy part of the /9-spectrum of '^H P^J2B[] . In the following we will perform all 
calculations in scheme A, but all bounds on CP-violating observables which we will derive 
hold in both schemes. 

Short-baseline neutrino oscillation experiments are sensitive to Am^ 
with a distance L between neutrino source and detector such that 



Ar/igi L 



< 1 



and 



Am43 L 



< 1 



(3.9) 



2p 2p 

As discussed in the previous section, with these assumptions on the neutrino mass spectrum 
there are no effects of CP violation in SBL neutrino oscillations. On the other hand, long- 
baseline neutrino oscillation experiments are planned to be sensitive to the "atmospheric 
neutrino range" 10~^eV^ < ^^Ij ~ 10~^eV^. In scheme A, the probabilities of Ua 
and Ua — ^ 1^/3 transitions in LBL experiments are given by 



p(LBL,A) 



Upi Ki + Up2 K2 exp 



. Am2i L 
' 2p 



k=3A 



(3.10) 



P,- 



(LBL,A) 



U^l Ual 



Ua2 exp 



Amni L 

>. 

2p 



k=3A 



(3.11) 



respectively. Matter effects are not included in these formulas which have been obtained 
from Eqs.( p.2|) and (|2.3|) , respectively, by taking into account that in LBL experiments 
Aml^L/2p ^ 1 and averaging out the oscillating terms with phases much larger than 27r 
{AmljL/2p > 2tt for A; = 3, 4 and j = 1,2). 

From Eqs.(|3J0|), ( pi]) and (^ it follows that the CP-odd asymmetries in LBL 

experiments in the schemes A and B are given by 



r)(LBL,A) 



n(LBL,B) 



r(A) 

r(B) 
'a/3 



sm 



sm 



2p 



Aml^L 



2p 



(3.12) 
(3.13) 
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with the oscillation amplitudes (see Eg. ( |2.14| )) 



/i^^=/a/3;12 and /^^ = /a/3;34 • (3.14) 

In the following we will only study scheme A and drop the superscript A (Z^^"* = for 
the reasons mentioned above in the context of Eq. ( p.8|) . 

Finally, we want to mention that the phases of the products of elements of the mixing 
matrix whose imaginary parts give the quantities Jg^, Jgr and /^^ are not independent |p8|| , 
as can be seen from the obvious relation 

arg [U,, U:, U;, U^2] = arg [U,, U;, Ul, t/.a] - arg [U,, U;, Ul, U,^] . (3.15) 

which is valid if Je^, ler and 1^^ are all different from zero. Hence, a measurement of le^, 
ler and I^r can give information on only two independent linear combinations of the three 
CP-violating phases which are possible in the four-neutrino schemes. In order to obtain 
information on the values of all the three CP- violating phases it is necessary to measure also 
some of the other Iai3-k/s. 



IV. CP VIOLATION IN THE SCHEMES WITH FOUR NEUTRINOS 



In this section we discuss bounds on the vacuum quantities ( p.l4|) . For the reasons 



explained in Sect. ^ we confine ourselves to scheme A in the derivations of the bounds. As 
shown in Appendix the unitarity of the mixing matrix implies the "unitarity bound" 

\Ia/3\ < fiCa,Cf3) (4.1) 

where f{x,y) is the continuous function 

( fi=xy for 2{1 - x){l - y) > xy 

f{x,y) = { (4.2) 
i /2 = 2[(x + y- 1)(1 - x)(l - y)]i/^ for 2(1 - x)(l - y) < xy 

defined on the unit square 0<X<1,0<?/<1. In Fig. [1| we have drawn a contour plot 
of the function f{x,y), which is helpful for the determination of the maximal allowed value 
for f{ca,C/3) when and/or are bounded. The dotted line in Fig. ^ is the borderline 
g{x) = 2(1 — x)/(2 — x) between the regions where f = fi and / = f2- Note that / is 
continuous along this borderline. 

In order to determine the maxima of f{x,y), the following considerations are useful 
(for the details consult Appendix 0). Increasing x at fixed y, the function / increases 
monotonously from / = at x = 0, until the straight line yi{x) = 2 — 2x (1/2 < x < 1) 
depicted in Fig. |I] is reached. There, the value of / is given by / = yy/1 — y. After 
this intersection, the function / decreases monotonously to / = at x = 1. From the 
symmetry f{x,y) = f{y,x), it follows that for fixed x and increasing y the function / 
increases monotonously from / = 0at?/ = 0to/ = xa/1 — x when the straight line 
1/2(2;) = 1 — x/2 (0 < X < 1) depicted in Fig. |l] is crossed. After this intersection, / 
decreases monotonously to f = at y = 1. 
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The absolute maximum of the function / (see Appendix |CD hes at the intersection of the 
hues yi and y2 and is given by /max = 2/3-\/3 ~ 0.385. Therefore, from the unitarity of the 
mixing matrix we have an absolute maximum for \Iai3\'- 

<^ ^ 0.385. (4.3) 

With the help of Fig. |l|, one can see that Eq.( [4.1|) with the constraints (|3.3| ) on Cg and 
implies that 

/2(«°,l/2(a°)) = a°(l-a°)^/^ for < > aO/2 , 

l^e.|<<; (4.4) 



The solid curve in Fig. shows the limit |/e^| < ctg 1^1 — a° with obtained from the 90% 
CL exclusion plot of the Bugey — > t'e experiment. The dash-dotted curve in Fig. |^ 
represents the improvement reached with the lower part of Eq. ([4.4|) at the values of Am^ 
for which a° < with a° obtained from the 90% CL exclusion plots of the CDHS |^ 



and CCFR [i3| t'^ — *■ experiments. From this figure one can see that the upper bound 



for |/e^| is very small (|/e^j| ^ 4 x 10 ^) for Am^ in the wide range (|2.18 



The bound represented by the solid curve in Fig. ^ is valid also for |JeT-|, because there 
is no experimental information on c^. 

For |/^r|; again by inspection of Fig. one can see that Eq.( |4.1j ) with the constraint 
( p.3|) on implies that 

I VI < /2(1 - aj,y2(l - <)) = (1 - aj) (4.5) 

The solid curve in Fig. |^ represents the corresponding bound obtained from the 90% CL 
exclusion curves of the CDHS and CCFR z/^ experiments. For Am^ < 0.3 eV^ 
there are no experimental data and therefore iVlmax ~ 0.385 by virtue of Eq.( [4.3|) . 



Taking into account the expression ( p.7| ) for Aa-js, in both schemes A and B we have also 
the "amplitude bound" (for the proof of this inequality, see Appendix 0) 



'a/3 1 



< ^\/Aa-t3{'^CaC^-A^.^). (4.6) 



With the upper bound A^.^ < we obtain 



1 



'a/3 1 



Ca Cf3 for A° > 2 Cc, . 

For |/e/i|, with the constraints (|3.3|), the inequality ( |4.7| ) becomes 



1 



for A0^,>2a0. 



(4.8) 



the Bugey pi 
E776 



The dashed curve in Fig. || shows the hmit ( [4.8|) obtained using the 90% exclusion plots of 
z/e experiment for the determination of 0° and the BNL E734 |45|, BNL 
experiments for the determination of ^J^.g- 



g and CCFR |3 



One can see 

that the upper bound for |Je^| is extremely small (|/e^i| ^ 10~^) for Am^ in the wide range 



( p.l8| ), which includes the LSND-allowed range ( |1.3| ). 

Since the constraints (|3.3|) do not put an upper bound on the possible values of and 
Cr, in the case of the inequality ( |4.7| ) becomes 



(4.9) 



The dashed curve in Fig. || shows the limit ( [4. 91) obtained using the 90% exclusion plots of 
the FNAL E531 ||8] and CCFR |59| experiments for the determination of 

The shadowed regions in Figs. H and ^ correspond to the range ( |1.3|) of Am^ allowed 
at 90% CL by the results of the LSND and all the other SBL experiments. From Fig. |^ it 
can be seen that, taking into account the LSND signal, |J^t-| could be close to the maximal 
value 2/3 allowed by the unitarity of the mixing matrix. 



V. THREE MASSIVE NEUTRINOS 



Though not all present indications in favour of neutrino mixing can be taken into account 
in scenarios with mixing of three massive neutrinos, it is nevertheless interesting to investi- 
gate also this case with the methods developed in this paper. Let us assume for definiteness 
that of the two differences of squares of neutrino masses one is relevant for SBL oscillations 
and the other one for LBL oscillations (see also Refs. p8|-p0|,^). These assumptions give 
rise to the following two three-neutrino mass spectra 



(I) 



LBL 



m\ < 1712 ^ ITT'S 
■' 



LBL 

and (II) mi <^ 777,2 < ^3 



(5.1) 



SBL SBL 

Furthermore, the results of the disappearance experiments allow to define the three regions 



(see Refs. 51,44,40 



(1) 
(2) 
(3) 



\Uek\^ > 1 



ek 



ek 



<a 



\U,k?<al 



(5.2) 



\Uij.k\ 



> 1 



with A; = 3 for the scheme I and k = 1 for the scheme II (for the definition of 0° and a°, 
see Eq.( p.5| )). The neutrino and antineutrino LBL oscillation probabilities in scheme I are 
given by 



p(LBL,I) 



)(LBL,I) 



Upi Ki + f//32 K2 exp 



2p 



Atti^ L 

f/Ji + f/j2 U^2 exp ( -I 



(5.3) 
(5.4) 
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Prom a comparison of Eqs. (|573|) and (|5.4|) with Eqs.( p.l(i| ) and ( |3.11|) , it is obvious that the 
CP-odd asymmetries D^^^^'^^ are given by the same formulas (|3.12| ) and ( p. 141) as in the 



4-neutrino case with superscript I instead of A. The transition probabihties in scheme II 
are obtained from the expressions ( |5.3|) and (|5.4|) by a cychc permutation of the indices: 
1 ^ 2 ^ 3 ^ 1. Therefore, as in the case of the schemes A and B for four neutrinos, the 
bounds on the CP-odd parameters Jq,^ are the same in the three neutrino schemes I and II. 

The methods for the derivation of the bounds on lajj which are described in the ap- 
pendices for the four-neutrino schemes ( p.l| ) are valid also in the case of mixing of three 
neutrinos. Obviously, the derivations of the four-neutrino case A (B) are carried over to the 
three-neutrino case I (II) if we put Uai = {Uai = and change the indices 2, 3, 4 — 1, 2, 3) 
for all a = e, fi, t. This implies that the amplitude bound (|A10|) applies also to the three 
neutrino schemes ( |5.1| ). In order to derive the unitarity bound from Eq.( |A10D , one must 



notice that from the unitarity of the 3x3 mixing matrix we have Aa;p = 4(1 — Ca)(l — cp) 
and Ca + Cf3 > 1. Hence, the equality sign applies in Eq . (p^) and there is no such distinction 
as defined by Eq. ([B^) . Consequently, by simple substitution of Aa-p = 4(1 — Cq,)(1 — c^) in 
Eq.( Al(]| ), in the case of mixing of three neutrinos one obtains the unitarity bound ( p6|) , i.e. 



\Ia/3\ < f2{Ca,Cf3) . (5.5) 

Notice that f2{caiCp) vanishes on the unitarity boundary Ca + cp = 1. Since the maxima 
of f{ca,cp) and f2{ca,C[^) coincide and are reached for Cq = = 2/3 (see Appendix y), 
the absolute maximum |/a/3|max = 2/3^/3 of the 4-neutrino case extends its validity to three 
neutrino^. 

Purthermore, since in the 3-neutrino case the CP-odd asymmetries in different oscillation 
channels are connected by Eq. (|2.11| ), we have 

-^e/^ — Ifj.T — Ire • (5-6) 



In the following we will give LBL bounds for each of the regions ( p.2|) , along the lines of 
the previous 4-neutrino section. 

Region 1. With respect to SBL and LBL neutrino oscillations, the 3-neutrino schemes 
I and II in Region 1 correspond to the 4-neutrino schemes A and B, respectively, with the 
same bounds on | Je^| (Eqs. (|4^ ), ( [4.8| ) and Pig. 0). Because of Eq.(|5.6|) the stringent bounds 
on \Iefi\ given in Pig. ^ are valid also for |/^r|- 

Region 2. Actually, this Region is disfavoured by the results of the LSND experiment 
(see Refs. [pTj^^^ ) . Also the results of the atmospheric neutrino experiments [^-[1^ 



taken together with the results of the CHOOZ experiment indicate that this Region is 



disfavoured. Indeed, in this region is small (c^ < 0° + a°) and the atmospheric neutrino 
anomaly can be explained only with dominant ^ Ve oscillations, which are forbidden 
by the results of the CHOOZ experiment. Since these evidences could disappear when the 



This value corresponds to the maximal value of the Jarlskog parameter J [^ 37] for CP viola- 
tion in the Kobayashi-Maskawa matrix, | J|max = l/6\/3, with an additional factor 4 due to the 



definition ( p. 141 ). 
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results of future more accurate experiments will be available, let us discuss the bounds on CP 
violation with the methods described in the appendices. They are given by the amplitude 
bound 



1 



< 1, V^°;e (4 - (5.7) 



and the unitarity bound 



Ie,\ < /2(1 - al 1 - aj) = 2 Jalal (l - - . (5.8) 



Both bounds are less restrictive than the corresponding bounds in Region 1. By Eg. ( p. 61) 
these bounds hold in all neutrino transition channels. 

Region 3. In this region, where Ce > 1 — and < aj^, the full set of atmospheric neu- 
trino data cannot be explained in the framework discussed here |26[. Applying nevertheless 



our methods for obtaining bounds on CP violation, we get the results 



and 




\U < \ \JaI,^ (4<-Ay (5.9) 



il-alY'^ for a0></2, , , 



2 [(a°-a°)(l-aj)a°]^/^ for a" < aO/2 . 



The amplitude bound is more stringent than the one in Region 2, but less restrictive than the 
one in Region 1 and in the 4-neutrino schemes A and B (for < aj^). From Eq.( |5.6| ) it follows 
that the bounds (5^) and ( 5.10 ) are valid also for the parameter |J^^| that characterizes the 



CP-odd asymmetry in the Ur channel. 

Summarizing the results of this section, we have shown that in regions 2 and 3 the strong 
bounds of region 1 and of the 4-neutrino schemes A and B on |/e^| are somewhat relaxed. 
However, regions 2 and 3 are disfavoured by present hints for neutrino oscillations. 

Let us also emphasize that the solar neutrino problem cannot be explained by neutrino 
oscillations in the three-neutrino schemes considered in this section. Hence, we regard them 
as remote possibilities. 



VI. MATTER EFFECTS AND CP VIOLATION 

Since in LBL neutrino oscillation experiments the neutrino beam travels a long distance 
through the earth's crust, matter effects influence the neutrino oscillation probabilities. The 
effective Hamiltonians in the flavour basis of neutrinos and antineutrinos in the case of 
mixing of four neutrinos are given, respectively, byQ 



Since active and sterile neutrinos are present in the schemes under consideration, both charged- 
current and neutral-current interactions contribute to the effective Hamiltonians. 
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H. = ^ [UM^U^ + diag (acc, 0, 0, a^c)) , (6.1) 

Hu = ^ (U*M^U^ - diag (ace, 0, 0, unc)) , (6.2) 

where acc and a^c are given by 

acc = 2V2G^iVep ~ 2.3 x lO^^eV^ ^ ^ ( 

\3gcm~'^y V 



.IGeV 



(6.3) 



ttNc = V2GFNnP ^^acc , (6.4) 

respectively, and = diag(m^, 777-2, m|, 777.4). Here is the Fermi constant, Ne and Nn 
are the electron and neutron number density, respectively, and p is the density of matter. 
With an average density of approximately 3 g cm~^ in the lithosphere we get 



2, -"-^^^i" iTtalj- 

Therefore, large matter effects are to be expected for baselines L > 1000 km. 

In the following we apply the simplifying approximation of constant electron and neutron 
number densities, which is rather accurate in the case of LBL experiments. In order to obtain 
the neutrino oscillation probabilities in matter we have to diagonalize the Hamiltonians ( |6.1|) 
and ( |6.2| ) with unitary matrices U' and U' , respectively, leading to the eigenvalues ej/2p of 

and ej/2p of Hp. Thus, we have 

H, = U'—U'^ and Hp = U'*—U'^ (6.6) 
2p 2p 

with the diagonal matrices e = diag(ei, €2, 63, 64) and e = diag(ei, 62, es, q) for neutrinos and 
antineutrinos, respectively (in the limit acc,o-NC ^ of vanishing matter effects we get 
U',U' U and ej,ej m]). 

Analogously to the definition of laisjk in Eq.( |2.14j ), it is useful to define 



j;^^,., ^ 4 Im [f/V f/- U':, U'p,] and 1'^^.^^, ^ 4 Im [^V f/^ U'^,] . (6.7) 

The usefulness of these definitions lies in the fact that for the proof of the existence of CP 
violation due to neutrino mixing the quantities I'a/3-jk ^'^^ -^ais-jk '^^ g'ooc? as la/s-jk- In other 
words, Iai3;jk = for all values of the indices a, k if and only if I'aipi-jik> = i^'a'/s'-j'k' — 0) 
for all values of a', k' . Let us prove this statement for /^/Jyfc (an analogous proof holds 
in the case of I'ap-jk)- Fii'st we assume that lap-jk = for all values of a,P,j,k. This is 
possible only if the mixing matrix U can be written a.s U = e'^'^Ue^'^, where p and a are real 
diagonal matrices and U is real. Since the matter part in the Hamiltonian ( |6.1| ) is real and 
diagonal, the effective Hamiltonian ( |6.1| ) can be written as 

= e*"^ (UM^W + diag(acc, 0, 0, a^c)) e"^^ = d'H^e-'" , (6.8) 
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where Hy is a real and symmetric matrix. Hence, the matrix U' can be written as U' = e*''f/' 
where U' is reaL This form of U' imphes that /^^.^^ = for all values of a,(3,j, k. In order 
to prove the inverse statement, i.e. that Ia/3;jk = for all values of k if I'^ipi-jik' = 

for all values of a', j', /c', we note that 

f/M^f/t = U'eU'^ - diag(acc, 0, 0, a^c) ■ (6.9) 

Now the roles of U and U' are exchanged. Hence, the same reasoning as before leads to the 
inverse statement. 

The possibility of finding evidence for CP violation in the lepton sector through the 
quantities I'afi;jk ^^"^ -^a/Syfc leads us to the search for methods that could allow to extract 
these quantities from the transition probabilities of neutrinos and antineutrinos measured 
in long-baseline oscillation experiments. 



The formula in Eq. (|2.12|) for the probability of ^Va^^^^p transitions is adapted to the 



matter case by the substitutions U ^ U' and I ^ I' for neutrinos and U ^ U' and I ^ I' 
for antineutrinos: 

Py.^y, = T.\uL,m,/ + 2j2MK,u';,^^^^^^^ (6.10) 

j k>j ^ k>j ^ 

p..^., = E + 2 E MKp';p':kU',k] ¥ - ^ E ^ > (e-n) 

j k>j ^ k>j ^ 

with ekj = Cfc — ej and ekj = — ej. From these two equations it is clear that in matter the 
transition probabilities of neutrinos and antineutrinos are different even if CP is conserved, 
i.e. if all the quantities I'ap-jk ^ap-jk equal to zero. Hence, simple measurements of the 



asymmetries (|2.8|) do not allow to obtain direct information on CP violation Pq-|30|. This 
is due to the fact that the matter contribution to the effective neutrino and antineutrino 
Hamiltonians ( |6.1|) , ( |6.2|) is not CP-symmetric. However, since the matter contribution to 
the effective neutrino (antineutrino) Hamiltonian is real and the matter density is symmetric 
along the path of the neutrino beam in terrestrial long-baseline experiments, matter effects 



are T-symmetric . In other words, there is no difference in the matter contributions to the 



'i/q— and '^/3— ^-'i^a channels and a difference of the corresponding transition probabilities 
can only be due to a fundamental violation of T in the lepton sector. Since the CPT theorem 
implies that a violation of T is equivalent to a violation of CP, we are lead to explore the 
possibility to obtain direct evidence of CP and T violation in the lepton sector through 
measurements of the T-odd asymmetries 

Ta/3 = Pua-^U/s ~ Pyp^UoL ^^'^ = Ppa^Up ~ Pu/S^Ua (6.12) 

in long-baseline oscillation experiments. In the case of a constant matter density along 
the neutrino path (which is a good approximation for terrestrial LBL experiments with a 
baseline shorter than about 4000 km), from Eqs.( |6.10| ) and ( |6.11D it is straightforward to 



obtain the following expressions for the T-odd asymmetries: 

Ta(3 = E ^Ujk sin ^ and = E ^afs-jk sin ^ . (6.13) 

k>j ^ k>j 
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It is clear that finding T^/j and/or Tap different from zero would be a direct evidence of T 
(and CP) violation in the lepton sector independent from matter effects. 

The T-odd asymmetries ( |6.12| ) cannot be measured in the accelerator LBL experiments 
of the first generation [pl| -|23|| because the initial beam will contain almost exclusively 



In order to have some possibility to measure the T-odd asymmetries it will be necessary to 
wait for the second generation of accelerator LBL experiments, as those that will use + z/g 
and + Ve neutrino beams from a muon collider ||^, 



Let us now consider the four- neutrino schemes ( |3.1j ), whose phenomenology of CP vi- 
olation in long-baseline neutrino oscillation experiments in vacuum has been discussed in 
Sections |TT| and |V[ We will discuss explicitly only the neutrino T-odd asymmetries T^^ in 
the scheme A, but all the conclusions are valid also for the antineutrino T-odd asymmetries 
TajB in scheme A and for the neutrino and antineutrino T-odd asymmetries in scheme B 
(with the exchange of indices given in Eq. (|3.8| )). 

In Ref. we have shown that apart from small corrections of order acc I ■, can 
decompose U' into 

U' = UR with R=[\'^ ^ \ (6.14) 
where -Ratm and -Rsun are 2x2 unitary matrices. The block structure of R implies that 

i=l,2 j=l,2 

and therefore 

Ca^Yl l^"^l'= E (6.16) 

i=l,2 j=l,2 

Since the bounds on Jq^ = Ia/3;i2 derived in Section |I^ depend only on quantities of the type 
( |6.15| ) and ( |6.16| ), we arrive at the interesting conclusion that the upper bounds on Iaf3 are 



also valid for 

^'afS = ^'af3;12 ^ud /^^ = 1^13-12 ■ (6-17) 

In this sense the upper bounds represented by the curves in Figs. and ^ include matter 
effects. 

Using a method analogous to that employed for the derivation of Eq.( p.l(]| ), for the 
probability of z/q, — > vp LBL transitions in matter we obtain 

= P'piKi + u',,u':,e-^*\' + \u',,u':, + %,u':,e-^-\' 



1 1 

+ 2^a/3si^'^+ 2 ^"Z^^^^'^' (6.18) 



with the definitions 
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= , ^ = -^"Z? = ^"/3;34 • (6-19) 

The expressions for the quantities 62 — ei and 64 — £3 in terms of the neutrino mixing 
parameters and of the matter density have been derived in Ref. . 



From Eq.( |6.18| ), for the neutrino T-odd asymmetries we obtain the expression 

Tai3 = I'afj sin (f) + J'^p sin uj . (6.20) 

Therefore, in matter the T-odd asymmetry T^^ depends not only on the parameter J^^ 
relative to the atmospheric sector (see Eq. ( |6.14| ) and the definitions (|6l^ ) and ( |6.17| )), but 
also on the parameter J^^ relative to the solar sector. In Ref. |Q we have shown that in 



the case of accelerator LBL experiments the maximal value of the parameter 00 is given by 

3 (i(j(jL 3 



2 2p 2 



V2Gi.Ar,L = 8.6 X 10"^ ( I (6.21) 

\ 1km / 



for p = 3gcm^^. Notice that the value of cUmax does not depend on the neutrino energy, but 
only on the propagation distance L. From Eq.( |6.21| ) one can see that the contribution of 
the term in Eq. (|6.2Cl|) proportional to sinu could be relevant for L > 100 km and therefore 
cannot be neglected in the analysis of the results of LBL experiments. 

As we have seen above, the upper bounds on 1^/3 = Iai3;i2 are also valid for J^^ and /^^. 
An analogous reasoning leads to the conclusion that the upper bounds on the parameters 

Jaf3 = IaP;U (6.22) 

are also valid for J^^ and J^^ = Hence we are lead to the investigation of the upper 

bounds for the parameters Jap- Following the methods presented in the Appendices 0-0, 
one can derive the unitarity bound 



\Jap\<f{l-c^A-cp) (6.23) 



and the amplitude bound 



1 



\Jap\<- \jA^,fi [4 (1 - c«) (1 - ca) - A^,fi\ . (6.24) 



Taking into account the constraints (|3.3|) on Ce and obtained from the exclusion curves 
of SBL reactor and accelerator disappearance experiments, the unitarity bound for | Je^| is 
given by 

, , ( f2{al,y2{al)) = al{l-alY^' for a° > aO/2 , 

\Jeii\ < \ (6.25) 
l/2«l-a°) = 2[(aj-a°)(l-aj)a0]^/^ a" < aO/2 . 

The numerical value of this bound, obtained from the 90% CL exclusion plot of the Bugey 
||1|] z/e — > z/e experiment and from the 90% CL exclusion plots of the CDHS and CCFR 



^ z/^ — >■ z/^ experiments, is shown by the solid curve in Fig. ^. The dashed curve in Fig. 



H represents the amplitude bound 
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for <2al, 
for AO^, > 2 a" , 



(6.26) 



obtained from the 90% CL exclusion plots of the CDHS H and CCFR |3 



experiments and from the 90% CL exclusion plots of the BNL E734 BNL E776 



and CCFR ^z/^^j— >-Ve experiments. The shadowed region corresponds to the range ( p..3|) 
of Am? allowed at 90% CL by the results of the LSND and all the other SBL experiments. 
From Fig. |^ it can be seen that I Je;.| < 10"^ for Am^ > 0.27 eV^ which is an interval that 
includes the LSND-allowed range (p.. 31). 



For \Jur\ we have the unitarity bound 



\J,r\<f2{al,y^{al))=al^l-al 



(6.27) 



This limit is more stringent than the corresponding one for |/^^| given in Eq.( [4.5|) . Its 



numerical value obtained from the 90% CL exclusion plots of the CDHS [E2| and CCFR H3 



experiments is shown by the solid curve in Fig. ^. The dashed curve in Fig. |^ 



represents the value of the amplitude bound 



for Al^^ < 2 a" , 
for Al^,>2al, 



(6.28) 



obtained from the 90% CL exclusion plots of the CDHS H and CCFR || 



experiments and from the 90% exclusion plots of the FNAL E531 and CCFR |5S 
experiments. From Fig. || one can see that |J^^| < 0.25 for Am? > 0.3 eV , 
including the LSND-allowed range ([L^ ), and \ J^r\ ^ 8 x 10~^ for Am^ > 0.5 eV^. 

Since 1 — Ce is large and there is no constraint on the value of (and the available 
information on Ae^r is rather poor), the parameter {Jerl is only subject to the unitarity 
bound 



|Je.| </2(l- a°, 2/2(1 -al^)) = (l-a°) 



(6.29) 



This bound is much less stringent than the corresponding one for {lerl, represented by the 
upper function in Eq.( [4.4| ) and depicted as the solid curve in Fig. ^ Since < 4 x 10~^ for 
Am^ in the wide range (|2.18|) , we obtain the upper bound iJerl ^ 0.2, which is about half 
of the unitarity limit 2/3^3 ^ 0.385. 

From the bounds depicted in Figs. H and ^, it is clear that the observation of a non-zero 
T-odd asymmetry T^g (and 7)je) in future LBL experiments is a very difficult task. On the 
other hand, from Figs. |^ and ^ one can see that the T-odd asymmetry T^^ (and T^r) could 
be rather large, close to the maximal value allowed by the unitarity of the neutrino mixing 
matrix. Also the T-odd asymmetry Ter (and Ter) could be rather large, but not more than 
half of the maximal value allowed by the unitarity of the neutrino mixing matrix and only 
if the matter effect is important and enhances the contribution of Jg^. Hence, we conclude 
that long-baseline ^ (and z/^ ^ Ur) experiments are favoured for the investigation of 
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CP (and T) violation in the lepton sector if the four-neutrino schemes are reahzed in 
nature. 

We would also like to mention that the CP-odd parameters are relevant for the 
CP-odd asymmetries that could be measured by extremely long-baseline (ELBL) oscillation 
experiments with neutrino beams propagating in vacuum, for which Am\iL/2p ^ 27r and 
A'm1^L/2p ~ 1 (in scheme A), 



P,(ELBL) J ■ Am^^L 

Da-p = Jap sm — , (6.30) 

Zp 



although we do not know if it will ever be possible to make such experiments. 

Concluding this Section, we briefly discuss the matter effects in the case of the 3-neutrino 
schemes considered in Section ^ We consider, for simplicity, only scheme I, but analogous 
conclusions are valid in scheme II. The 3-neutrino scheme I corresponds to the 4-neutrino 
scheme A, with the difference that the solar sector of the mixing matrix is absent, i.e. 
Ua4: = for a = e, /i, r and i?sun = 1 (apart from negligible corrections of order acc/Afn"^)- 
Hence, in the 3-neutrino scheme I the neutrino T-odd asymmetries T^/j (see Eg. ( p.l2| )) are 
given by Eq.( |6.20|) with J^^ = 0. Using the same reasoning as that employed in the case of 



the 4-neutrino schemes, one can see that the upper bounds on I^jj derived in Section |V| in 
the case of the 3-neutrino schemes are also valid for /^^ and /^^. Hence, T^^ is very small if 
the neutrino mixing parameters lie in region 1 and is less suppressed in regions 2 and 3. 



VII. CONCLUSIONS 

In this paper we have considered possibilities to reveal effects of CP-violation in the 
lepton sector in future accelerator long-baseline (LBL) experiments (K2K ||21]], MINOS p2 
ICARUS m and others 



At present there are three experimental indications in favour of neutrino oscillations 
which correspond to three different scales of neutrino mass-squared differences: the solar 
neutrino deficit, the atmospheric neutrino anomaly and the result of the LSND experiment. 
These indications and the negative results of numerous short-baseline (SBL) neutrino ex- 
periments can be accommodated in the two four-neutrino schemes A and B presented in 



Working in the framework of the neutrino mixing schemes A and B, we have derived 
constraints on the parameters Jq,/? (a,/3 = e,/i, r) (see Eq. (|3.14|) ) that characterize CP 



violation in ^i^a^^i^l LBL neutrino oscillation experiments in vacuum. These parameters are 
given in terms of the quantities la/B-jk which appear in the general CP-odd asymmetries (see 
Eqs. (|2.13| ) and ( |2.14D ) as lap = Ia(3;i2 in scheme A and lap = Iap;M in scheme B. We have 
developed methods for deriving upper bounds on the parameters lap from the data of SBL 
experiments which can be applied not only to the schemes A and B but to arbitrary schemes 
with any number of neutrinos. We have shown that the CP-odd parameter Jg^ is bounded by 
\Iep\ < 10~^ (see Fig. 0) and a similar suppression applies to ler- On the other hand, sizable 
CP violation can be expected in oscillations. The CP-odd parameter relative to this 

channel could be close to its maximally allowed value |/^r|max = 2/3-\/3 ~ 0.385, resulting 
from the unitarity of the mixing matrix (see Fig. 13). 



18 



For LBL accelerator experiments the matter background is important. In this case the 
parameters la/sjk for neutrinos (antineutrinos) are replaced by the CP-violating parameters 
^ai3;jk i^'ap;jk) which includc matter effects (see Eg. (|6.7D). Using a physically motivated 
approximate method of incorporating matter effects |5^, we have demonstrated that the 
quantities |/^^| are bounded by the same functions of the SBL mass-squared difference 

Am? as \Iai3\ apart from terms of order acc I ■ Therefore, the bounds depicted in Figs. |^ 
and ^ apply also to the corresponding /^^ and /^^. However, although finding /^^ 7^ and/or 
7^ would prove the existence of CP violation in the neutrino mixing matrix f/, the 
knowledge of the parameters /^^ and/or 7^^ cannot easily be transformed into information on 
Jq,/3 and thus f/, because both sets of parameters are related in complicated, non-linear way 
involving also the "matter potentials" ace-, o-nc and the mass-squared difference relevant 
for LBL neutrino oscillations. 

We have also shown that in matter a second CP-violating parameter J^^ (given by 
J'^p = 1^^.34 in scheme A and J^^ = /^^.i2 scheme B) appears in the z/q, — >• transition 
probability. Its contribution is only significant if the oscillation phase uj defined in Eq. ( |6.19|) 
is sufficiently large. An evaluation of the maximal value that the phase uj can assume 
in accelerator LBL experiments shows that the contribution of J^^ could be relevant for 
baselines longer than ~ 100km (see Eq.( |6.2lD ). We have argued that the parameter J'^^ 
(and the analogous parameter J^^ for antineutrinos) is subject to the same bounds as Ja^ 
(with Ja/3 = /q,/3;34 in scheme A and Ja/s = Ia/3;i2 in scheme B). These bounds are presented 
in Eqs.( |6.25| )- (|6.28|) and their numerical values obtained from the results of disappearance 
and appearance SBL neutrino oscillation experiments are shown by the curves in Figs. ^ 
and ^ There is no analogue of the parameters J^^ and J^^ in the 3- neutrino case. 

Since a measurement of the CP-odd asymmetries Dajs defined in Eq.( p.8|) does not 
allow to obtain direct information on CP violation if matter effects are important, we have 
considered the long-baseline T-odd asymmetries Ta;/3 (Ta;/?) defined in Eq.( |6.12|) . Since the 
matter contribution to the effective neutrino and antineutrino Hamiltonians ( |6.1| ), ( |6.2| ) is 
real and the matter density is symmetric along the path of the neutrino beam in terrestrial 
LBL experiments (to a good approximation it is even constant for baselines shorter than 
~ 4000 km) the matter effects are T-symmetric. Therefore, the T-odd asymmetries Tq,.^ 
(Tq,.^) are only different from zero if CP is violated in the lepton mixing matrix U. We 
have shown that in the four-neutrino schemes A and B the T-odd asymmetries depend on 
the parameters J^^ {I'a/3) ^.nd J^^ iJ'a/s) (^^^ Eq. (|6.20|) ). Hence, they are subject to the 
constraints derived from the results of SBL experiments. 

In conclusion, we have shown that in the four-neutrino schemes A and B the channels 

'i^^^'z/e are disfavoured for the search of CP- violation effects in future LBL oscillation exper- 
iments, the channels 'z/e^=>'^v could allow to reveal relatively large CP- violating effects where 
matter plays an important role and the channels *z/^^=^*^v could show CP-violating effects as 
large as is allowed by the unitarity of the neutrino mixing matrix. 
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APPENDIX A: DERIVATION OF THE AMPLITUDE BOUND 



In this appendix we discuss the derivation of the "amphtude bound". The starting point 
is the quantity 

hp = A\m[U^^U;,U*^^Up2\ (a 7^/9), (Al) 

which determines the CP-odd asymmetry in the four-neutrino scheme A. The same bound 
can be derived in the scheme B with the change of indices (|3.8|) . 

It is obvious that Jq,^ is invariant under the phase transformations 

U^j ^ e^^^ f/,,- , ^ e'^^ Up, , (A2) 

where the jj are arbitrary phases. Thus the elements Uaj can be taken to be real. Taking 
into account the definition f|3.2|), we can write 



Uaj = y/c^eY^ with j = l,2 (A3) 



,(1) 

and the orthonormal basis 

e(i) (6) = (cos e, sin 6) , e^^) {6) = {- sin 6, cos 6) . (A4) 
We expand t/gj with respect to this basis as 

p=l,2 

where pi and p2 are complex coefficients such that 

Y.\P,\' = 1. (A6) 

p=l,2 

With the help of Eqs.(|X|)-(^ we easily find 

Iai3 = 2caCj3 siu 29 Im(p^p2) = 2 Cq, \pi I a/1 — sin 29 sin 6 , ( A7) 

where 6 is the phase of plp2- 

The parameter \pi\ is connected to the oscillation amplitude A^-p and the parameters 
Ca, cp. In fact, from Eqs. ( |A3|) and (|A5|) we have 

2 



A. 



a;l3 



i=i,2 



4c„c/3|pip, (A8) 



which implies = \J Aa-p/ AcaCp. Inserting this expression in Eq. ([AT]) , we obtain 



lap = 2 V ^"^'f^ ~ "^"'^^ ^^"^ 



and thus we arrive at the "amplitude bound" 

|/a/3| < A^,p{Ac^cp- A^,p). (AlO) 

Let us stress that this derivation is based only on the obvious inequality 

|sin2^sin5| < 1. (All) 
Since Cq,, C/3 and Aa-^p do not restrict 9 and 5, the bound ( [AlOp is the optimal one. 



20 



APPENDIX B: DERIVATION OF THE UNITARITY BOUND 



Up to now we did not use the unitarity of the mixing matrix. Taking this fact into 
account will allow us to obtain an upper bound on \Ia/3\ depending solely on Cq, and cp. 
The unitarity of the mixing matrix tells us that 



j=l,-2 i=3,4 



(Bl) 



This relation allows to write the oscillation amplitude Aa-j3 in the two forms of Eq.( |3.7] ). 
Using the Cauchy-Schwarz inequality, one can see that (in the scheme A) 



i=3,4 



< 



\i=3,4 / \j=3,4 / 



4(l-c„)(l-c^). (B2) 



The right-hand side of the inequality (|A10|) , as a function of Aa-f3, reaches its maximum, 
CaCj3, at (here we do not take into account possible experimental information on Aq,.^) 

{Aa-i3)o = 2CaCf3 . (B3) 

Consequently, if the condition 

2(l-c,)(l-c^) >CaC^ (B4) 
is satisfied, the upper bound ( |B^ ) on Aq,.^ is larger then (v4q,.^)o. In this case we have 

1 4/3 1 < CaCf}. (B5) 

If the condition ( P^ ) is not fulfilled, the upper bound ( P^ is smaller than {Aa-i3)o and has 
to be inserted for Aa-fs into Eq. ( |A10| ), leading to 



\Iaf3\ <2^{Ca + Cp-l){l- C„) (1 - C/3) . 

Thus, we arrive at the "unitarity bound" 

with the function 

fi=xy for 2{1 - x){l - y) > xy , 

h = 2[{x + y-l){l-x){l- y)Y'^ for 2(1 - x){l - y) < xy , 
defined on the unit square 0<x<l,0<|/<l. The function 

. . 2(l-a;) 



f{x,y) 



(B6) 



(B7) 



[B8) 



9[x) 



2-x 



(B9) 



represents the borderline separating the two regions in the definition of the function ( |B8D . 
It is clear from our derivation (and also easy to check) that / is continuous along this 
borderline. 
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APPENDIX C: DISCUSSION OF THE FUNCTION / 



Since we do not have definite experimental values of Cq, and C/j, but only bounds on these 
quantities (see Eqs.( |3.3|) and ( |3.4| )), which define allowed rectangles in the square < Cq, < 1, 
< C/3 < 1, we are interested in the behaviour of / in order to evaluate the unitarity bound. 

From the partial derivative of / in the region y > g{x), 



one can see that, at fixed y, the function / increases monotonously from x = to the point 
X = 1 — y/2, where the partial derivative in Eq. (|Cl|) is zero. The points x = 1 — ?//2 lie on 
the straight line yi{x) = 2 — 2x. In the range l-~y/2<x<l the function / decreases 
monotonously. Taking into account the symmetry f{x,y) = f{y,x), we see that at fixed x 
the function / increases monotonously from y = Oto the point y = l — x/2, where the partial 
derivative of / with respect to y is zero. These points lie on the straight line y2{x) = l — x/2. 
Beyond this line / decreases monotonously. Note that both straight lines lie in the range of 



Figure |I] shows a contour plot of the function f{x,y), together with the lines yi and y2 
which intersect at the point 




(x{2-2x-y) 



(CI) 



/2. 



2 



(C2) 




At this point both partial derivatives of / are equal to zero and therefore the point ( |C2|) 
corresponds to the absolute maximum of /, given by 




(C3) 



This number constitutes the absolute upper bound for \Iaf3\- 



22 



REFERENCES 



[1] J.H. Christenson, J.W. Cronin, V.L. Fitch and R. Turlay, Phys. Rev. Lett. 13, 138 
(1964). 

[2] W. Bernreuther, CP Violation in Weak Decays and Elsewhere^ talk given at the 20th 
Johns Hopkins Workshop, Heidelberg, June 1996, preprint [hep-ph/9701357| ; Andrzej 
J. Buras, CKM Matrix: Present and Future, talk given at the Symposium on Heavy 
Flavours, Santa Barbara, July 7- July 11, 1997, preprint [hep-ph/9711217| . 

[3] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). 

[4] B.T. Cleveland et al, Nucl. Phys. B (Proc. Suppl.) 38, 47 (1995). 

[5] K.S. Hirata et al, Phys. Rev. D 44, 2241 (1991). 

[6] GALLEX Coll., W. Hampel et al, Phys. Lett. B 388, 384 (1996). 

[7] J.N. Abdurashitov et al, Phys. Rev. Lett. 77, 4708 (1996). 

[8] K. Inoue, Talk presented at TAUP97, September 7-11, 1997, Laboratori Nazionali del 
Gran Sasso, Assergi (Italy). 

[9] Y. Fukuda et al, Phys. Lett. B 335, 237 (1994). 
[10] R. Becker-Szendy et al, Nucl. Phys. B (Proc. Suppl.) 38, 331 (1995). 
[11] W.W.M. Allison et al, Phys. Lett. B 391, 491 (1997). 

[12] K. Martens, Talk presented at the International Europhysics Conference on High Energy 
Physics, 19-26 August 1997, Jerusalem, Israel ( [http:/ /www.cern.ch/hep97/abstract/-| 
tpalO.html) ; E. Kearns, Talk presented at the Conference on Solar Neutrinos: News 



About SNUs, 2-6 December 1997, Santa Barbara, California (|http:/ /doug-pc.itp.ucsb." 
I edu/online/snu/l) . 

[13] C. Athanassopoulos et al., Phys. Rev. Lett. 77, 3082 (1996). 

[14] GALLEX Coll., P. Anselmann et al, Phys. Lett. B 285, 390 (1992); P.I. Krastev and 

S.T. Petcov, ibid. 299, 99 (1993); N. Hata and P.G. Langacker, Phys. Rev. D 50, 632 

(1994); G. Fiorentini et al, ibid. 49, 6298 (1994). 
[15] V. Barger, R.J.N. Phillips and K. Whisnant, Phys. Rev. Lett. 69, 3135 (1992); P.I. 

Krastev and S.T. Petcov, Phys. Rev. Lett. 72, 1960 (1994). 
[16] S.P. Mikheyev and A.Yu. Smirnov, Yad. Fiz. 42, 1441 (1985) [Sov. J. Nucl. Phys. 42, 

913 (1985)]; II Nuovo Cimento C 9, 17 (1986); L. Wolfenstein, Phys. Rev. D 17, 2369 

(1978); ibid. 20, 2634 (1979). 
[17] S.M. Bilenky and B. Pontecorvo, Phys. Rep. 41, 225 (1978). 
[18] S.M. Bilenky and S.T. Petcov, Rev. Mod. Phys. 59, 671 (1987). 

[19] C.W. Kim and A. Pevsner, Neutrinos in Physics and Astrophysics, Contemporary 
Concepts in Physics, Vol. 8, (Harwood Academic Press, Chur, Switzerland, 1993). 

[20] V. Barger et al, Phys. Rev. D 22, 1636 (1980); S.M. Bilenky and F. Niedermayer, Yad. 
Fiz. 34, 1091 (1981) [Sov. J. Nucl. Phys. 34, 606 (1981)]. 

[21] Y. Suzuki, Proc. of Neutrino 96, Helsinki, June 1996, edited by K. Enqvist et. al., p. 
237 (World Scientific, Singapore, 1997). 

[22] MINOS Coll., D. Ayres et al, NUMI-L-63, February 1995. 

[23] ICARUS Coll., P. Cennini et al, LNGS-94/99-I, May 1994. 

[24] NOE Coll., G.C. Barbarino et al, INFN/AE-96/11, May 1996; AQUA-RICH Coll., 
T. Ypsilantis et al, LPC/96-01, CERN-LAA/96-13; OPERA Coll., H. Shibuya et al, 
CERN-SPSC-97-24. 



23 



[25] S.M. Bilenky, C. Giunti and W. Grimus, Proc. of Neutrino 96, Helsinki, June 1996, 
edited by K. Enqvist et. al., p. 174 (World Scientific, Singapore, 1997). 

[26] S.M. Bilenky, C. Giunti and W. Grimus, preprint |hep-ph / 9607372] , to appear in Z. Phys. 
C. 

[27] N. Okada and O. Yasuda, Int. J. Mod. Phys. A 12, 3669 (1997). 
[28] M. Tanimoto, Phys. Rev. D 55, 322 (1997). 

[29] J. Arafune and J. Sato, Phys. Rev. D 55, 1653 (1997); J. Arafune, M. Koike and J. 

Sato, Phys. Rev. D 56, 3093 (1997). 
[30] H. Minakata and H. Nunokawa, preprint TMUP-HEL-9704 ( |hep-ph/9703^ ; preprint 

TMUP-HEL-9705 ( |hep-ph/ 97062811) . 
[31] T.K. Kuo and J. Pantaleone, Phys. Lett. B 198, 406 (1995). 
[32] S. Geer, preprint Fermilab-PUB-97/389 ( |hep-ph/9712'290|) . 

[33] R.N. Mohapatra, Talk presented at the Workshop on Physics with the First Muon 
Collider, Fermilab, 6-9 November 1997, preprint liep-ph/9711444| . 

[34] M. Doi et al, Phys. Lett. B 102, 323 (1981); L. Wolfenstein, Phys. Lett. B 107, 77 
(1981); S.M. Bilenky et al, Nucl. Phys. B 247, 61 (1984); B. Kayser, Phys. Rev. D 30, 
1023 (1984). 

[35] S.M. Bilenky, J. Hosek and S.T. Petcov, Phys. Lett. B 94, 495 (1980); M. Doi et al, 

Phys. Lett. B 102, 323 (1981). 
[36] C. Jarlskog, Phys. Rev. Lett. 55, 1039 (1985); Z. Phys. C 29, 491 (1985); C. Jarlskog 

and R. Stora, Phys. Lett. B 208, 268 (1988); C. Jarlskog, Proc. of CP Violation, Ed. C. 

Jarlskog, Advanced Series in High Energy Physics Vol. 3 (World Scientific, Singapore, 

1989), p. 3. 

[37] D.D. Wu, Phys. Rev. D 33, 860 (1986); I. Dunietz, O.W. Greenberg and D.D. Wu, 

Phys. Rev. Lett. 55, 2935 (1985). 
[38] J.D. Bjorken and I. Dunietz, Phys. Rev. D 36, 2109 (1987); I. Dunietz, Ann. Phys. 184, 

350 (1988). 

J.T. Pehoniemi and J.W.F. Valle, Nucl. Phys. B 406, 409 (1993); D.O. Caldwell and 
R.N. Mohapatra, Phys. Rev. D 48, 3259 (1993); Z. Berezhiani and R.N. Mohapatra, 
Phys. Rev. D 52, 6607 (1995); J.R. Primack et al, Phys. Rev. Lett. 74, 2160 (1995); E. 
Ma and P. Roy, Phys. Rev. D 52, R4780 (1995); R. Foot and R.R. Volkas, Phys. Rev. D 
52, 6595 (1995); E.J. Chun et al, Phys. Lett. B 357, 608 (1995); J.J. Gomez-Cadenas 
and M.C. Gonzalez-Garcia, Z. Phys. C 71, 443 (1996); S. Goswami, Phys. Rev. D 55, 
2931 (1997); A.Yu. Smirnov and M. Tanimoto, Phys. Rev. D 55, 1665 (1997); E. Ma, 
Mod. Phys. Lett. A 11, 1893 (1996). 



[40 
[41 
[42 
[43^ 
[44^ 
[45^ 
[46 
[47 
[48^ 
[49^ 



S.M. Bilenky, C. Giunti, C.W. Kim and S.T. Petcov, Phys. Rev. D 54, 4432 (1996). 

B. Achkar et al, Nucl. Phys. B 434, 503 (1995). 

F. Dydak et al, Phys. Lett. B 134, 281 (1984). 

I.E. Stockdale et al, Phys. Rev. Lett. 52, 1384 (1984). 

S.M. Bilenky, A. Bottino, C. Giunti and C.W. Kim, Phys. Rev. D 54, 1881 (1996). 

L.A. Ahrens et al, Phys. Rev. D 36, 702 (1987). 

L. Borodovsky et al, Phys. Rev. Lett. 68, 274 (1992). 

A. Romosan et al, Phys. Rev. Lett. 78, 2912 (1997). 

N. Ushida, Phys. Rev. Lett. 57, 2897 (1986). 

K.S. McFarland et al, Phys. Rev. Lett. 75, 3993 (1995). 



24 



[50] S.M. Bilenky, C. Giunti and W. Grimus, preprint |hep-ph/ 97102091 , Phys. Rev. D 57 
(1998) in press. 

[51] S.M. Bilenky, A. Bottino, C. Giunti and C.W. Kim, Phys. Lett. B 356, 273 (1995). 
[52] M. ApoUonio et al, preprint |hep-ex/ 971 1002 . 



25 



FIGURES 



FIG.|l[ Contour plot of the function f{x,y) given in Eq. (|4.2|) . The dotted hne is the 
borderhne g{x) = 2(1 — x)/(2 — x) between the regions where / = /i and / = f2- The two 
sohd hues represent the functions yi{x) = 2 — 2x and y2{x) = 1 — x/2. 



FIG.^ Upper bound for the parameter |Je^| which characterizes the CP-odd asymmetry 
in the z/^ channel for the SBL parameter Am^ in the range 10-1 eV^ < Am^ < lO^eVl 
The sohd curve represents the upper function in Eq. ([4.4|) and is obtained from the 90% CL 
exclusion plot of the Bugey experiment. The dash-dotted curve improves the 

solid curve where aj^ < flg/^ (the lower function in Eq. (|4.4|) ). It is obtained from the 90% 
CL exclusion plots of the Bugey Ve — ^ experiment and the CDHS and CCFR z/^ 
experiments. The dashed curve is obtained from the 90% CL exclusion plots of the Bugey 
Ue — »■ t'e experiment and the BNL E734, BNL E776 and CCFR z/^ — > and — > z/g 
experiments (see the upper function in Eq.( ^4.8| )). The shadowed region corresponds to the 
range ( |1.3| ) of Am^ allowed at 90% CL by the results of the LSND experiment. The solid 
curve represents also an upper bound for |/er|. 



FIG.^ Upper bound for the parameter |/^t-| which characterizes the CP-odd asymmetry 
in the z/^ Vr channel. The solid curve is obtained from the 90% CL exclusion plots of 
the CDHS and CCFR — > z/^ experiments (see Eq. ([4.5|) ). The dashed curve is obtained 
from the 90% CL exclusion plots of the FNAL E531 and CCFR experiments (see 

Eq.( [4.9|) ). The shadowed region corresponds to the range (|1.3|) of Am^ allowed at 90% CL 
by the results of the LSND experiment. 



FIG.|[ Upper bound for the parameter | Je^| (see Eq.( |6.22| )). The solid curve represents 
the unitarity bound (|6.25|) and is obtained from the 90% CL exclusion plots of the Bugey 
z/g — > z/e experiment and those of the CDHS and CCFR z/^ — > z/^ experiments. The dashed 
curve represents the value of the amplitude bound ( |6.26| ) obtained from the 90% CL exclusion 
plots of the CDHS and CCFR z/^ experiments and those of the BNL E734, BNL E776 
and CCFR z/^ v^, and z/^ z/g experiments. The shadowed region corresponds to the 
range (|1.3|) of Am^ allowed at 90% CL by the results of the LSND experiment. 



FIG.^ Upper bound for the parameter | J^,-! (see Eq. (|6.22| )). The solid curve represents 
the unitarity bound ( |6.27D obtained from the 90% CL exclusion plots of the CDHS and CCFR 
experiments. The dashed curve depicts the value of the amplitude bound (|6.28|) 
obtained from the 90% CL exclusion plots of the CDHS and CCFR z/^ — z/^ experiments and 
those of the FNAL E531 and CCFR z/^ Uj. experiments. The shadowed region corresponds 
to the range (|l]3|) of Am^ allowed at 90% CL by the results of the LSND experiment. 
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Figure [J 
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